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Abstract 



Within the framework of the Top-color-assisted Technicolor (TC2) model, we 
calculate the new physics contributions to the branching ratios B(B S ^ — > 77) 
and CP violating asymmetries rQ P (B s ^ — > 77) in the QCD factorization 
based on the heavy-quark limit mj S> Aqcd- Using the considered parameter 
space, we find that (a) for both B s — ► 77 and B& — > 77 decays, the new 
physics contribution can provide a factor of two to six enhancement to their 
branching ratios, (b) for the B s — > 77 decay, its direct CP violation is very 
small in both the SM and TC2 model, and (c) the CP violating asymmetry 
r~^ p {Bd — > 77) is around the ten percent level in both the SM and TC2 model, 
but the sign of CP asymmetry in the TC2 model is different from that in the 
SM. 
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I. INTRODUCTION 



As is well known, the rare radiative decays of B mesons induced by the quark decay 
b — > 57 (q — d, s) are very sensitive to the flavor structure of the standard model (SM) and 
to new physics beyond the SM. Both inclusive and exclusive processes, such as the decays 
B — > X s 7, X s 77 and B s ^ — > 77, have been studied in great detail [1-11]. 

The inclusive decay B — > X s 7 is measured experimentally with increasing accuracy [12]. 
The world average as given by the 2002 Particle Data Group ( PDG2002 ) [13] is 

B(B -> X sl ) = (3.3 ± 0.40) x 10~ 4 , (1) 

which is quite consistent with the next-to-leading order (NLO) standard model prediction 
[4] 

B(B -> X sl ) TH = (3.29 ± 0.34) x 10" 4 , (2) 

Obviously, there is only small room left for new physics effects in flavor- changing neutral 
current processes based on the b — > s transition. In other words, the excellent agreement 
between SM theory and experimental data results in a strong constraint on many new physics 
models beyond the SM. 

Within the SM, the electroweak contributions to b — > S77 and B — > 77 decays have been 
calculated some time ago [1]; the leading-order QCD corrections and the long-distance con- 
tributions were evaluated recently by several groups [6,14]. The new physics corrections were 
also considered, for example, in the two-Higgs doublet model [15,16] and the supersymmetric 
model [17]. 

On the experimental side, only upper limits (90%C.L.) on the branching ratios of B s ^ — > 
77 are currently available 

B(B S -> 77) < 1.48 x 10~ 4 , [18] (3) 
£(£ d ^ 77) < 1.7 x 10" 6 , [19] (4) 

which are roughly two orders above the SM predictions [1,6,8,9]. These radiative decays 
are indeed very interesting because (a) these decays have a very clean signal where two 
monochromatic energetic photons are produced precisely back-to-back in the rest frame of 
B meson; (b) these exclusive decays also allow us to study the CP violating effects as the 
two photon system can be in a CP-even or CP-odd state; (c) since B s — > 77 depends on 
the same set of Wilson coefficients as B — > X s 7, its sensitivity to new physics beyond the 
SM complements the corresponding sensitivity in B — > X s 7; and (d) the smallness of the 
branching ratios can be compensated by the very high statistics expected at the current B 
factory experiments and future hadron colliders. 

In this paper, we present our calculation of branching ratios and CP-violating asymme- 
tries for rare exclusive decays B s d — > 77 in the framework of the Top-color-assisted Techni- 
color (TC2) model [20] by employing the QCD factorization based on the heavy-quark limit 
m b > Aqcd [8,9]. 

This paper is organized as follows: In Sec. II, we give a brief review about the SM 
predictions for the branching ratios and CP asymmetries of B St d — > 77 decays. In Sec. Ill, 
we present the basic ingredients of the TC2 model, and evaluate the new penguin diagrams. 
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After studying the constraint on the TC2 model by considering the data of B° d - B° d mixing 
and B — > X s 7 decay, we find the Wilson coefficients Cj and Cg with the inclusion of the 
new physics (NP) contribution. In Sec. 4, we show the numerical results of branching ratios 
and CP-violating asymmetries for B s d — > 77 decays. The discussions and conclusions are 
included in the final section. 



II. B s ,d -»• 77 DECAYS IN THE SM 

In this section, based on currently available studies, we present the formulae for exclusive 
decay B s ^ — > 77 in the framework of SM. 

A. Effective Hamiltonian for inclusive b — > S77 decay 

We know that the quark level processes b — > sj, S77 and the exclusive decays B s ^ — > 77 
have a close relation. Up to the order the effective Hamiltonian for the decay b — > S77 

is identical to the one for 5 — > X s 7 transition [1,7] 

W e //(6 - *7) = We//(6 - ^77) + O . (5) 



m 



w , 



This can be understood by either applying the equation of motion [21] or by applying an 
extension of Low's low energy theorem [22]. 

Up to corrections of order l/rriyy, the effective Hamiltonian for b — > 577 is just the one 
for b — > S7 and takes the form 



ft = !f£ V A (s) 

V ^ p=u,c 



CMQI + CMQI+ ]T CiMQi 

i=3,...,8 



(6) 



where = V^* V^, * s the Cabibbo-Kobayashi-Maskawa (CKM) factor. And the current- 
current, QCD penguin, electromagnetic and chromomagnetic dipole operators are given by 



Qi = {s a Pp)v-A{Pi3b a )v-A, (7) 
Ql = (sp) V -A(pb)v-A, (8) 

Q3 = (sb) v -A J2 (qq)v-A, (9) 

q=u,d,s,c,b 

Qi = (S a bp)v-A J2 (QpQa)v-A, (10) 
q=u,d,s,c,b 

Q5 = (sb) v ^ A J2 (qq)v+A, (11) 

q=u,d,s,c,b 



1 For the numbering of operators Qi j2 > we use the convention of Buras et al. [3] throughout this 
paper. 
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where a and (3 are color indices, 



Qq = (s a b/3) V ~A ^2 (Q/3Qo)v+A, 
q=u,d,s,c,b 

g 

Qi = s a a' u/ (m b R + m s L) b a F^ v , 
Qs = 4^ s a a^ (m b R + m s L) T^bp G% , 



(12) 

(13) 
(14) 



a 



1. 



, 8 labels SU (3) c generators, e and g s refer to 
the electromagnetic and strong coupling constants, and L,R — (1 =p 7s)/2, while and 
G°^ v denote the photonic and gluonic field strength tensors, respectively. In Qi^ the terms 
proportional to m s are usually neglected because of the strong suppression m 2 /m 2 . The 
effective Hamiltonian for b — > c/77 is obtained from Eqs.(6) - (14) by the replacement s — > d 
The Wilson coefficients Cj(/i) in Eq.(6) are known currently at next-to- leading order 
(NLO) [2,3]. Within the SM and at scale m\y, the Wilson coefficients Ci{m w ) at the 
leading order (LO) approximation have been given for example in [3], 



Ci(m w ) = 
C 2 {m w ) = 1, 

C 7 {m w ) = 



1,3,4,5,6), 



-7x t 



5x 2 t 



8x 3 



2x1 



24(1 - x t f 



C 8 (m w ) 



—2x t — 5x 2 + x 



8(1 -x t f 



4(1 -x t y 

3x 2 t 

4(1 -x t y 



■log[x t ], 



(15) 
(16) 

(17) 
(18) 



where x t = mf/m^. 

By using QCD renormalization group equations [3], it is straightforward to run Wilson 
coefficients Ci(m w ) from the scale \x = 0(m w ) down to the lower scale \x = 0(m b ). The 
leading order results for the Wilson coefficients Ci(fi) with ji ps m fe are of the form [3] 



C,-(A*)=EM a< (J = !,-••, 6), 

i=i 

C 7 (/i) = i^C^) + I (?/ 4 / 23 - r/ 6 / 23 ) C 8 (m w ) + £ ^ , 



i=i 



i=l 



where 77 = a s (mw)/oi s (fj), and the magic numbers are [3] 

(14/23, 16/23, 6/23, -12/23, 0.4086, -0.4230, -0.8994, 0.1456) , 



(2.2996, -1.0880, -3/7, -1/14, -0.6494, -0.0380, -0.0185, -0.0057) , 

(0.8623, 0, 0, 0, -0.9135, 0.0873, -0.0571, 0.0209) , 

/0 1/2 -1/2 \ 

1/2 1/2 

-1/14 1/6 0.0510 -0.1403 -0.0113 0.0054 

-1/14 -1/6 0.0984 0.1214 0.0156 0.0026 

-0.0397 0.0117 -0.0025 0.0304 

\0 0.0335 0.0239 -0.0462 -0.0112/ 



(19) 
(20) 
(21) 



(22) 
(23) 
(24) 



(25) 
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The numerical results of the LO Wilson coefficients Ci(fi) obtained by using the input 
parameters as given in Table I are listed in Table II for (j, = m b /2, m b and 2m b , respectively. 

B. B Sj d — ► 77 decays in the SM 

Based on the effective Hamiltonian for the quark level process b — > 5(^)77, one can write 
down the amplitude for B s d — > 77 and calculate the branching ratios and CP violating 
asymmetries once a method is derived for computing the hadronic matrix elements. There 
exist so far two major approaches for the theoretical treatments of exclusive decay B — > 77. 

The first approach was proposed ten years ago and has been employed by many authors 
[1,6]. Under this approach, one simply evaluates the hadronic element of the amplitudes 
for one-particle reducible (1PR) and one-particle irreducible (1PI) diagrams, relying on a 
phenomenological model. One can work, for example, in the weak binding approximation 
and assume that both the b and the light q quarks are at rest in the B q meson [23]. From 
the heavy quark effective theory(HQET), for instance, one can also assume that the velocity 
of the b quark coincides with the velocity of the B q meson up to a residual momentum of 
h-QCD- Both pictures are compatible up to corrections of order (AgcWmf,) [23]. One typical 
numerical result obtained by employing this approach is 

B(B S -> 77) « (2 - 8) x KT 7 (26) 

after inclusion of LO QCD corrections [23]. There are also many works concerning the 
estimation of the long distance contributions to B — > 77 decay [14]. 

In the first approach, one has to employ hadronic models to describe the B q (q = s, d) 
meson bound state dynamics. It is thus impossible for one to separate clearly the short- 
and long-distance dynamics and to make distinctions between the model-dependent and 
model-independent features. 

The second approach was proposed recently by Bosch and Buchalla [8,9]. They analyzed 
the B s ^ — > 77 decays in QCD factorization approach based on the heavy quark limit m b ^> 
Aqcd- Under this approach, one can systematically separate perturbatively calculable hard 
scattering kernels from the nonperturbative B-meson wave function. Power counting in 
AqcD/nT'b allows one to identify leading and subleading contributions to B — > 77. In this 
paper, we will employ the Bosch and Buchalla (BB) approach to calculate the Technicolor 
corrections to B s ^ — > 77 decays. 

From Refs. [8,9], one knows that (a) only one 1PR diagram [ Fig. 1(a) ] contributes at 
leading power; (b) the most important subleading contributions induced by the 1PR [Fig. 
1(b) ] and 1PI diagrams [Fig. 1(c) ] can also be calculated; and (c) the direct CP violation 
of Bd — > 77 can reach the 10% level. 

The amplitude for the B — > 77 decay has the general structure [8] 

A(B -> 7 (*i,ci)7(*2,C2)) = ^/4^IW F - iA -^l°>- (27) 

Here F^ v and F pv = e fJ,uXp F Xp /2 are the photon field strength tensor and its dual with 
£ oi23 _ _^ rj^g Draric hing ratio of B q — > 77 decay with q = s, d is then given by 
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B[B q -> 77) = tb, 



G 2 F ml 



fl 



288tt 3 



I A 



(28) 



where (7p is the Fermi constant, a em is the fine structure constant, tb q is the lifetime of B q 
meson, and and are the mass and decay constant of the B q meson, respectively. 
The values of all input parameters are listed in Table I. 

The matrix elements of the operators Qi in Eq.(6) can be written as 



(7(ei)7(e2)|Qi|S) = f B CdtTnt) $b(0^2, 

jo 



(29) 



where the are the polarization 4-vectors of the photons, = $51 is the leading twist 
light-cone distribution amplitude of the B meson, and T^(£) is the hard-scattering kernel 
describing the hard-spectator contribution. 

By explicit calculations as were done in Ref. [8], the quantities A± in Eq.(28) are of the 
form 



A + = X^A U + + X^A C + , 
A_ = X^Al + X^AI, 



(30) 
(31) 



with 



A P : 



-C 7 ^ + (C 5 + 3C 6 ) 



1 . . 1 

j 9 (l)-3 



Ab o 



2g(z c ) + -g(l) 



(C 4 — Cq 



\g{z c ) + ^(1) 



where z p = m^/ml for p = u,c, and 



20^ A „ 16 „ 
+ —C 3 + 4C 4 - —C 5 , 



(32) 



(33) 



g{z) = -2 + 4z 



Li-. 



1 - y/l-Az + i 



+ Li 2 



it , 



1 + ^/l-Az + i 



1€ , 



(34) 



and Li 2 (x) is the dilogarithm function. It is easy to see that A\ = A c + , but A u _ ^ A c _. The 
function g(z) has an imaginary part for < z < 1/4, while g(0) = —2 andg(l) = 2(rr 2 — 9)/9. 

The first term of A± is the leading power contribution from the 1PR diagram [ Fig. 1(a) ] of 
the penguin operator Q 7 , the remaining terms of A v ± represent the subleading contributions 
from the 1PR diagram [ Fig. 1(b) ] with the operator Q 7 where the second photon is emitted 
from the b quark line, and from the 1PI diagram [ Fig. 1(c) ] induced by insertion of four- 
quark operators Qi. From the formulas as given in Eq.(28) and Eqs.(30) - (33), we find the 
numerical results of the branching ratios in SM 



B(B S 
B(B d 



77) 
77) 



1 2 +2A 



(AA B )^(A/i) ± 0.3(A/ B J ± 0.02(A 7 ) 



x 10" 



3.2^g(AA B )lg : g(A^)lS;g(A/ fl(l )lSi(A7) 



x 10" 



(35) 
(36) 



where the central values of branching ratios are obtained by using the central values of input 
parameters as given in Table I, and the errors correspond to AX B = ±0.15 GeV, m b /2 < 
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fj, < 2mt), AfB d = Afs s = ±0.03 GeV, respectively. For the CKM angle 7, we consider 
the range of 7 = (60 ± 20)° according to the global fit result [13]. Obviously, the dominant 
errors are induced by the uncertainty of hadronic parameter A#, the renormalization scale 
/j, and decay constant f Bq . The error induced by A7 is about 30% for B d decay, but very 
small for B s decay. The errors due to the uncertainty of other input parameters are indeed 
very small and can be neglected. 

Now we consider the CP violating asymmetries of B S)d — > 77 decays. Following the 
definitions of Ref. [8], the subscripts ± on A± for B — > 77 decay denote the CP properties 
of the corresponding two-photon final states, while A± refer to the CP conjugated amplitudes 
for the decay B — > 77 (decaying b antiquark). Then the deviation of the ratios 



\A± 


2 _ 


A± 


2 


\A± 


2 + 


A ± 


to | 



(37) 



from zero is a measure of direct CP violation. Since A p + = A+, Tq P is always zero. For r^p 
of Bd — > 77 decay, however, it can be rather large. By using the central values of input 
parameters as given in Table I and assuming 7 = 60°, we find 



It is easy to see that the direct CP violating asymmetry for B s — > 77 decay is small, 
~ 1%, and cannot be detected by experiments. For B d — > 77 decay, however, its CP 
violation can be rather large, around —10% for 7 ~ 60°. But the much smaller branching 
ratio is a great challenge for the current and future experiments. 

In Fig. 2, we show the CKM angle 7— and /i-dependence of r^ P {B d — > 77). The 
dotted, short-dashed and solid curves show the SM predictions of r^ P {B d — > 77) for /i = 
mb/2,mb and 2m&, respectively. The CP violating asymmetry even can reach —17% for 
CKM angle 7 50°, the value preferred by the global fit [24] and by the analysis based on 
the measurements of branching ratios of B — > Kn decays [25]. The value of r^, p (B d — > 77) 
here is the same as that given in Ref. [8] for ji = nib, but opposite with what was given in 
Ref. [8] for ji = nib/ 2 and 2nib, respectively. 



In this section, we calculate the loop corrections to B Sjd — > 77 decays in TC2 model. 



Apart from some differences in group structure and/or particle contents, all TC2 mod- 
els [20,26] have the following common features: (a) strong Top-color interactions, broken 
near 1 TeV, induce a large top condensate and all but a few GeV of the top quark mass, 
but contribute little to electroweak symmetry breaking; (b) technicolor [27] interactions are 
responsible for electroweak symmetry breaking, and extended technicolor (ETC) [28] inter- 
actions generate the masses of all quarks and leptons, except that of the top quarks; (c) 




(38) 
(39) 



III. B s ,d 77 DECAYS IN TC2 MODEL 



A. TC2 model 
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there exist top pions fr^ and 7r° with a decay constant Fq ps (40 — 50) GeV. In this paper 
we will chose the well-motivated and most frequently studied TC2 model proposed by Hill 
[20] to calculate the contributions to the rare exclusive B decays from the relatively light 
charged pseudo-scalars. It is straightforward to extend the studies in this paper to other 
TC2 models. 

In the TC2 model [20], after integrating out the heavy coloron and Z' ', the effective 
four-fermion interactions have the form [29] 

Ceff = Mj { ( K + Iff) ^ R ^ L + { K ~ ^) WhM^} , (40) 

where k = (gl/An) cot 2 9 and K\ = (gf/Air) cot 2 9', and My is the mass of coloron V a and 
Z 7 . The effective interactions of Eq. (40) can be written in terms of two auxiliary scalar 
doublets 0i and 2 . Their couplings to quarks are given by [30] 

C eff = \l4> L <l)ltR + \2^ L 4>^>R, (41) 

where A 2 = Att(k + 2«i/27) and A 2 = 4ti(k — K\/27). At energies below the top-color scale 
A ~ 1 TeV the auxiliary fields acquire kinetic terms, becoming physical degrees of freedom. 
The properly renormalized <f>i and 2 doublets take the form 



where 7r ± and 7r° are the top pions, H ±,a and A are the b pions, h t is the top Higgs boson, 
and Fq 50 GeV is the top pion decay constant. 

From Eq. (41), the couplings of top pions to t and b quark can be written as [20] 



TYv — 

i ttTT° + i t R b L 7T + + i —t L b R n + + H.c. 



(43) 



where m* t = (1 — e)m t and ml m I GeV denote the masses of top and bottom quarks 
generated by top-color interactions. 

For the mass of top pions, the current 1 — a lower mass bound from the Tevatron data 
is m^ > 150 GeV [26], while the theoretical expectation is m^ ~ (150 — 300GeV) [20]. For 
the mass of b pions, the current theoretical estimation (100 - 350)GeV 

and mfj = m 2 ^ + 2mf [30]. For the technipions nf and 7r^, the theoretical estimations 
are m ni > 50GeV and m na ps 200GeV [31,32]. The effective Yukawa couplings of ordinary 
technipions nf and 7r^ to fermion pairs, as well as the gauge couplings of unit-charged 
scalars to gauge bosons 7, Z° and gluon are basically model-independent, can be found in 
Refs. [31-33]. 

At low energy, potentially large flavor- changing neutral currents (FCNC) arise when the 
quark fields are rotated from their weak eigenbasis to their mass eigenbasis, realized by the 
matrices Ul,r for the up-type quarks, and by D L R for the down-type quarks. When we 
make the replacements, for example, 

b L ^D b L d d L + D b L s s L + D b L b b L , (44) 
b R - D^dn + D%s R + D bb b R , (45) 
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the FCNC interactions will be induced. In TC2 model, the corresponding flavor changing 
effective Yukawa couplings are 



i n + (D b L s i R s L + D b fi R d L ) + % H + (D b °t L s R + D%t L d R ) + H.c. . (46) 



For the mixing matrices in the TC2 model, authors usually use the "square-root ansatz" : 
to take the square root of the standard model CKM matrix (Vckm = U^D L ) as an indication 
of the size of realistic mixings. It should be denoted that the square root ansatz must be 
modified because of the strong constraint from the data of B° — B° mixing [30,34,35]. In 
the TC2 model, the neutral scalars H° and A can induce a contribution to the B® — 
(q = d, s) mass difference [29,30] 



AMr 7 m 



M Bq - UF^mir^ 



' < B Bq Fl, (47) 



where M Bq is the mass of B q meson, F Bq is the S g -meson decay constant, B Bq is the renor- 
malization group invariant parameter, and 5b q ~ \D b ^D b ^\. For the B d meson, using the ex- 
perimental measurement of AM Bd = (3.22±0.05) x lQ- l0 MeV [13] and setting F Q = 45GeV, 



B Bd F Bi = 200MeV, one has the bound 5 bd < 0.82 x 10~ 7 for m m < 600GeV. This is an 

important and strong bound on the product of mixing elements D 1 ^^. As pointed in [29], if 
one naively uses the square-root ansatz for both D L and D R , this bound is violated by about 
2 orders of magnitude. By taking into account above experimental constraint, we naturally 
set that D 1 ^ = for % ^ j. Under this assumption, only the charged technipions nf, 7T^ and 
the charged top pions ^ contribute to the decays studied here through penguin diagrams. 



B. Constraint on TC2 model from B — > X s ~f decay 

The constraint on both Dl and D R from the experimental data of B — > X s ^ decay is 
much weaker than that from the B° — B° mixings [29]. On the other hand, one can draw 
strong constraint on the mass of top-pion m# from the well measured B — > X s 7 decay by 
setting D b L d = V td /2, D b L s = V ts /2, F Q = 45 GeV and e = 0.05 ± 0.03. 

In this subsection, we firstly calculate the new physics contributions to the Wilson coeffi- 
cients Cj(mw) and C$(mw)- And then we draw the constraint on the mass m% by comparing 
the theoretical prediction of B(B — > X s 7) with the measured value as given in Eq.(l). 

The new photonic- and gluonic-penguin diagrams can be obtained from the corresponding 
penguin diagrams in the SM by replacing the internal W ± lines with the unit-charged scalar 
(ni ,7rf and 7r ± ) lines, as shown in Fig. 3. For details of the analytical calculations, one 
can see Ref. [36]. 

By evaluating the new 7-penguin and g/won-penguin diagrams induced by the exchanges 
of three kinds of charged pseudoscalars (tt ± , nf , nf), we find that 
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where y t = m|/((l - e)m t ) 2 , rj t = ml 1 /(em t ) 2 , & = m 2 8 /(em t ) 2 , while the functions H(x), 
K{x) and L(x) are 

, 22 - 53x + 25x 2 3a; - 8x 2 + 4x 3 
* (X) = 36(1 -x)3 + 6(1 l0gM ' (50) 

5 - 19s + 20a: 2 x 2 - 2x 3 
K{X)= 12(1 -x)3 2TT^F l0g[a:] ' (51) 

LW= 12(1 -,)■ + 2(1^ '° gM - <52> 

It is easy to show that the charged top-pion strongly dominate the new physics contri- 
butions to the Wilson coefficients C 7 (mw) and Cg(mw), while the technipions play a minor 
rule only, less than 5% of the total NP correction. We therefore fix the masses of nf and ir~^ 
in the range of m n = 200 ± 100 GeV and = 400 ± 100 GeV, as listed in Table III. At 
the leading order, the charged- scalar s do not contribute to the remaining Wilson coefficients 

Qi — Q&- 

When the new physics contributions are taken into account, the Wilson coefficients 
Cj(mw) and Cs(mw) can be defined as the following, 

C 7 {mw) Tot = C 7 {m w ) SM + C 7 (m w ) TC2 , (53) 
C 8 (m w ) Tot = C 8 (m w ) SM + C s (m w ) TC2 , (54) 

where C 7 ^f have been given in Eqs.(17), (18). Explicit calculations show that the Wilson 
coefficients Cj$ 2 have the opposite sign with their SM counterparts, and therefore they will 
interfere destructively. The QCD running of C 7 ot from the energy scale m w to fx w m b is 
the same as the case of SM. 

Using the NLO formulas as presented in Ref. [4] for the B — > decay, we find the 
numerical results for the branching ratios B(B — > X s 7) in both the SM and the TC2 model, 
as illustrated in Fig. 4, where we use the central values of input parameters as given in Table I 
and Table III. The three curves correspond to fx — rrib/2 (short-dashed curve), fx = m b (solid 
curve) and fx = 2m b (dot-dashed curve), respectively. The band between two horizontal 
dotted-lines shows the SM prediction B(B -> X s ^) = (3.29 ±0.34) x 10~ 4 [4], while the band 
between two horizontal solid lines shows the data, 2.5 x 10~ 4 < BiB — > X s ^f) < 4.1 x 10~ 4 
at the 2a level [13]. 

From Fig. 4 and considering the errors induced by varying m ni , m ns and e in the ranges 
as shown in Table III, the constraint on the mass of charged top pion is 

m # = 200 ± 30GeV, (55) 

which is a rather strong constraint on mj. 



IV. NUMERICAL RESULTS IN TC2 MODEL 

In this section, we present the numerical results for the branching ratios and CP violating 
asymmetries of B s ^ — > 77 decays in the TC2 model. 
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A. Branching ratios B(B St d — ► 77) in TC2 model 

Based on the analysis in previous sections, it is straightforward to present the numerical 
results. Our choice of input parameters are summarized in Table I and Table III. Using the 
input parameters as given in Table I and Table III and assuming 7 = (60 ± 20)°, we find 
the numerical results of the branching ratios 

B(B S -> 77) = [2.8l^(AA B )l^(A / ,)^; 8 7(A/ Bs )l^(Am # )] x 10" 6 , (56) 
B(B d -> 77) = [8.2^ 7 2 °(AA B )^|(A / ,)^; 3 7 (A/ B Jl^(Am # )] x 1(T 8 , (57) 

where the major errors correspond to the uncertainties of A\ B = ±0.15 GeV, m b /2 < pi < 
2mft, Af Bq = ±0.03 GeV and Amj = 30 GeV, respectively. 

Figures 5(a) and 5(b) show the charged top-pion mass and /i-dependence of the decay 
rates B{B S ^ — > 77), respectively. In these figures, the lower three lines show the SM 
predictions for \i = mfe/2(dotted line), ji = m^solid line) and ji = 2mb (short-dashed line). 
Other three curves correspond to the theoretical predictions of TC2 model. The new physics 
enhancement on the branching ratios and their scale and mass dependence can be seen easily 
from the figure. 

From the numerical results as given in Eqs.(56), (57), it is easy to see that the largest 
error of the theoretical prediction comes from our ignorance of hadronic parameter X B . We 
show such \ B dependence of branching ratios in Fig. 6 explicitly. The dotted and short- 
dashed curves in Fig.6 show the SM predictions for \x = mb/2 and /i = mj,, respectively. The 
dot-dashed and solid curves show the TC2 model predictions for \x = mb/2 and \x = mj,, 
respectively. The decay branching ratios decrease quickly, as \ B getting large for both SM 
and TC2 model. 

In order to reduce the errors of theoretical predictions induced by the uncertainties of 
input parameters, we define the ratio R(B q — > 77) with q = d, s as follows 

**->-^£. 

Using the central values of input parameters, one finds numerically that 

R(B S -> 77) = 2.34 ± 0.10(AA B )l 1 ; 2 7 °(A /J )l^(Am # ), (59) 
R(B d - 77) = 2.56 ± 0.02(AA B )l?;18(A^)l^(Am jf ), (60) 

where the errors correspond to AA B = ±0.15 GeV, mb/2 < fi < 2mb and Amj = 30 GeV, 
respectively. The dependence on input parameters f B , m 3 B , Gp and a em cancelled in the 
ratio R. 

In Figs. 7(a) and 7(b), we show the /x, mj and \ B dependence of the ratio R explicitly. 
It is easy to see from Fig. 7(b) that the strong A B -dependence of the individual branching 
ratios is now greatly reduced in the ratio R, but the strong /i-dependence still remains large. 
Obviously, the new physics enhancements to both branching ratios can be as large as a 
factor of two to six within the reasonable parameter space. 
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B. Direct CP violation of B S)d — ► 77 in TC2 model 



Now we calculate the new physics correction on the CP violating asymmetries of B sd — > 
77 decays. By using the input parameters as given in Tables I and III, we find the numerical 
results as follows: 



r- CP {B s -> 77) TC2 = -0.25lg;S(A^) ± 0.10(AA B ) ± 0.04(Am,)r^(A 7 ) x 10" 2 , (61) 



+0.10/ a ,,\ ±nin^\J-in rW a ™^+o.o7/ 



rap(S d - 77) iGi = [+6.5l&(A/i)3;i!(AA B ) ± l.l(Am*)l^(A 7 )j x 10~ 2 , (62) 

where the major errors are induced by the uncertainties of the corresponding input param- 
eters AA B = ±0.15 GeV, m b /2 < fi < 2m b , Amj = 30 GeV and A7 = ±20°, respectively. 

For the B s — > 77 decay, its direct CP violation is still very small after inclusion of new 
physics corrections. For the B d — > 77 decay, however, its CP violating asymmetry is around 
7% in TC2 model and depends on the hadronic parameter X B , the scale //, the CKM angle 
7 and the mass mj, as illustrated by Figs. 8 and 9. 

In Fig. 8 we draw the plots of the CP violating asymmetry r^ P (B d — > 77) versus the 
parameters fi, Xb and 7. The lower and upper three curves in Fig. 8 show the theoretical 
predictions of the SM and TC2 model, respectively. In Fig. 8(b), 7 = 60° is assumed. It 
is easy to see from Fig. 8 that the pattern of the CP violating asymmetry in TC2 model is 
very different from that in the SM. The sign of r' CP (B d — > 77) in TC2 model is opposite to 
that in the SM, while its size does not change a lot. Such difference can be detected when 
the statistics of the current and future B experiments becomes large enough. 



V. DISCUSSIONS AND SUMMARY 

In this paper, we calculate the new physics contributions to the branching ratios and CP- 
violating asymmetries of double radiative decays B Std — > 77 in the TC2 model by employing 
the QCD factorization approach. 

In Sec. II, based on currently available studies, we present the effective Hamiltonian for 
the inclusive B — > X s 7 and b — > 577 decays. For the evaluation of hadronic matrix elements 
for the exclusive B Syd — > 77 decays, we use Bosch and Buchalla approach to separate and 
calculate the leading and subleading power contributions to the exclusive decays under study 
from 1PR and 1PI Feynman diagrams. We reproduce the SM predictions for the branching 
ratios B(B sd — > 77) and direct CP asymmetries r^ p as given in Ref. [8]. 

For the new physics part, we firstly give a brief review about the basic structure of 
TC2 model, and evaluate analytically the strong and electroweak charged- scalar penguin 
diagrams in the quark level processes b — > s/d^ and b — > sg. We extract out the new 
physics contributions to the corresponding Wilson coefficients C 7 {m w ) and C s (mw)- Then 
we combine these new functions with their SM counterparts and run these Wilson coefficients 
from the scale \x = my/ down to the lower energy scale \x = 0{m b ) by using the QCD 
renormalization equations. From the data of B^ — B® mixing, we find the strong constraint 
on the "square-root ansatz". We also extract the strong constraint on the mass m% by 
comparing the theoretical predictions for the branching ratio B(B — > X s ^) at the NLO level 
with the experimental measurements. 
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In Sec. IV, we present the numerical results for B(B St d — > 77) and r CP (B St d — > 77) after 
the inclusion of new physics contributions in the TC2 model. 

1. For both B s — > 77 and Bd — > 77 decays, the new physics contribution can provide a 
factor of two to six enhancement to their branching ratios. The mj, fi and Xb depen- 
dences are also shown in Fig. 5. With an optimistic choice of the input parameters, 
the branching ratio B(B S — > 77) and BiBd —> 77) in the TC2 model can reach 1CT 5 
and 1CT 7 respectively, only one order away from the experimental limit as given in 
Eqs.(3), (4). With more integrated luminosity accumulated by BaBar and Belle Col- 
laborations, the upper bound on BiBd — > 77) will be further improved, and may reach 
the interesting region of TC2 prediction. 

2. For the B s — > 77 decay, its direct CP violation is very small in both the SM and TC2 
model. 

3. For the Bd — > 77 decay, however, its CP violating asymmetry is around ten percent 
level in both the SM and Tc2 model. But the pattern of CP violating asymmetry in 
TC2 model is very different from that in the SM, as illustrated in Fig. 8. 

As discussed in Ref. [37], the high luminosity option SuperBaBar suggests a total inte- 
grated luminosity of 10 ab -1 . For the branching ratio as given in Eq.(57), the number of 
observed Bd — > 77 events is then expected to be in the range of 50 — 150 in the TC2 model, 
and therefore measurable in the future. 
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TABLES 

TABLE I. Values of the input parameters used in the numerical calculations. All masses are 
in units of GeV. 



A 


A 




Gf 


O^em 


a s (M z ) 


0.847 


0.2205 


0.38 ± 0.08 


1.1664 x 10" 5 GeL- 2 


1/137.036 


0.118 


m w 


m t 


pole 
m b 


m pole 


m Bd 




80.42 


175 


4.80 ±0.15 


1.4 ±0.12 


5.279 


5.369 




I'b s 




11 MS 


r{B d ) 


r(B s ) 



: _ : : : m : : : : 

0.20 ±0.03 0.23 ±0.03 0.35 ± 0.15 0.225 1.542pa 1.461ps 



TABLE II. The LO Wilson coefficients Cj(//) in the SM obtained by using the central values 
of input parameters as listed in Table I. 



A* 


Ci 


c 2 


c 3 


c A 


c 5 


c 6 


c 7 


c 8 


m b /2 


-0.3500 


1.1630 


0.0164 


-0.0351 


0.0096 


-0.0467 


-0.3545 


-0.1649 


m b 


-0.2454 


1.1057 


0.0109 


-0.0254 


0.0073 


-0.0309 


-0.3141 


-0.1490 



2m b -0.1654 1.0664 0.0070 -0.0175 0.0052 -0.0200 -0.2801 -0.1353 



TABLE III. Values of the input parameters of TC2 model. All masses are in units of GeV. 







m% 


F n 




e 


200 ± 100 


400 ± 100 


200 ± 30 


120 


45 


0.05 ± 0.03 
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FIGURES 



s/d 



s/d 



' i 


•s/d 


(a) 


« i 


' b 



7 s/d 




(ft) 



FIG. 1. The leading power 1PR diagram (a) and subleading power 1PR diagram (b) of the mag- 
netic penguin operator and the subleading power 1PI diagram (c) of the four-quark operators 
Qi. The diagrams with interchanged photons are not shown. 
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FIG. 2. The CP violating asymmetry of (B^ — > 77) decay vs the CKM angle 7 and energy scale 
fi in the SM. The dotted, short-dashed and solid curves show the SM predictions for /i = mb/2, m\, 
and 2mfe, respectively. 
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9, 7 7 
FIG. 3. The typical photon- and gluon-penguin diagrams with W and charged-PGB exchanges 
(short-dashed lines) in the SM and TC2 models which contribute to B — > decays. The 

internal quarks are the upper type u, c and t quarks. 




FIG. 4. The branching ratios B{B — > X s j) in the SM and TC2 models as a function of mj. 
The band between two horizontal dashed-lines (solid lines) shows the SM prediction (world average 
of experimental measurements) as listed in Eqs.(l), (2). The short-dash, solid and dot-dash curves 
show the TC2 model predictions of the branching ratios for /j, = mb/2,mb and 2m&, respectively. 
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FIG. 5. Plots of branching ratios B(B S — ► 77) (a) and B{Bd — ► 77) (b) vs m^, setting Ab = 0.35 
and CKM angle 7 = 60°. The lower three lines in each diagram show the SM predictions for 
fi = mi,/ 2 (dotted line), (i = nif, (solid line) and [i = Im^ (short-dashed line). Upper three curves 
correspond to the theoretical predictions of TC2 model. 
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FIG. 6. Plots of branching ratios B(B S — ► 77) (a) and B(B& — ► 77) (b) vs A#, setting mj = 200 
GeV and CKM angle 7 = 60°. The dotted and short-dashed curves show the SM predictions for 
[i = rrib/2 and [i = nib, respectively. The dot-dashed and solid curves show the TC2 model 
predictions for (i = rrib/2 and fi = m^, respectively. 
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FIG. 7. The ratio of branching ratios R(B St d — ► 77) in the TC2 model. The three dotted and 
three solid curves show the ratios for B s — ► 77 and — ► 77 decays, respectively. In (b), we set 
= 170 GeV. 
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FIG. 8. The CP violating asymmetry of (B d — > 77) decay in the SM and TC2 model. The lower 
(upper) dotted, short-dashed and solid curves show the SM (TC2) predictions for \i = mb/2,mi } 
and 2mfe, respectively. 
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FIG. 9. The CP violating asymmetry of (B^ — > 77) decay vs mass and energy scale [i in 
TC2 model. The dotted, short-dashed and solid curves show the TC2 predictions for u, = m.b/2, vrib 
and 2rrib, respectively. 
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